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In the standard model (SM), the CP violation is introduced through a single phase in the CKM 
matrix. The neutral kaon system is one of the most precise channels to test how the SM theory 
describes the experiment data such as Ek accurately. The indirect CP violation is parametrized into 
Ek, which can be calculated directly using lattice QCD. In this calculation, the largest uncertainty 
comes from two sources: one is Bk and the other is V c b- We use the lattice results of Bk and 
exclusive V c b to calculate the theoretical estimate of Ek, which turns out to be 3 . 1 <7 away from its 
experimental value. Here, the error is evaluated using the standard error propagation method. 
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1. Introduction 

The neutral kaon system has two kinds of CP violation: indirect and direct CP violation. 
Indirect CP violation is parametrized into %. The experimental value of % is very well known [|l|] : 

e K = (2.228 ±0.011) x 1(T 3 xe*, e = 43.51 ±0.05. (1.1) 

We can also calculate Ek directly from the SM using Bk, V c b, and other input parameters, which 
are determined from other experiments and the SM theory. By comparing these two values, we can 
test a fundamental ansatz of the SM, the unitarity of the CKM matrix. 

Here, we calculate Ek directly from the SM using the known parameters with their errors in 
control. Two of the most important input parameters are Bk and V c b, which dominate the statistical 
and systematic uncertainty in Ek ■ During a past decade, lattice QCD has reduced the Bk error 
dramatically down to less than 5% level [||-||] as well as V c b [§]. 

There are two independent methods to determine V c b'. inclusive channels and exclusive chan- 
nel. There exists about 2 3a difference in V c b between inclusive and exclusive methods. We address 
this issue on how this have an effect on Ek- 

Here, we use the Wolfenstein parameters for the CKM matrix to calculate Ek mainly because 
they are convenient. 

Let us define a parameter A which test the unitarity ansatz directly as 

A = V-V^-I. 

Here, we calculate Ek and A to test the SM. The CKM matrix elements are accurate up to G(X 5 ) « 
3 x 10~ 4 . We use the standard error propagation method to estimate the errors. 



2. Review of the Neutral Kaon Mixing: e K 

The neutral kaon system forms a two dimensional Hilbert space. In this subspace, the time 
evolution of the neutral kaon state vectors can be described by the effective Hamiltonian // e ff 

i-\K(t))=H eS \K(t)), H efi - = M-i-. (2.1) 
at 2 

M (f ) is the dispersive (absorptive) part. Here, the dispersive part gives the mass eigenvalues of K$ 
and Ki, and the absorptive part represents the decay rates of K$ and K^. Let us take the basis with 
the CP even \K{) and odd \K 2 ) states, 

\K 1 ) = -L(\K°) + \t ) )), \K 2 ) = -L(\K°) -{?>)). (2.2) 

In this basis, the matrix elements can be written as the followings: 

M=( M \ +im \r=( r \ +i A. (2.3) 

y-im' m 2 \-iy r 2 J 

By construction of the formal perturbation theory in quantum field theory known as the Wigner- 
Weisskopf theory [Q], M and T are hermitian matrices. In addition, if we assume that the CPT 
invariance is exactly respected, then m' and / must be real. 
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Solving the eigenvalue problem with this matrix representation for H e ff, the eigenstates 

\K s ) = —l==(\K l ) + ~e\K 2 )), \K L )= J_ (\K 2 ) + e\Ki)) (2-4) 

V l + \ £ \ Vi + N 

have the mass M$,l and the decay rate Vs t, respectively. The small CP impurity £ satisfies the 
following equation 

S = S (0) (1 + S 2 ). (2.5) 

Here, the £(o) parameter is denned as 

£(0) = ~, ZJll =e M s i n 0(— j cot 0_^_) + ^(g3 ) (2.6) 

(0) (M 1 -M 2 )-i(r 1 -r 2 ) amk af k j ^ { 

where 

2AM K 

AM K = M L - M s , Ar K = r s -r L ,tand = —±. (2.7) 

AF K 

The solution of £ can be obtained by iteration. Since £( ) is of the order of 10~ 3 , we may expand e 
perturbatively as follows, 

£ = £(o) + £(o) + 2 ^fo) + 5£(o) H • (2-8) 



In Eq. (2.6), M\ 2 and T\ 2 can be safely replaced by the eigenvalues Ms,l and r$,L which are 



experimental observables. Note that this approximation makes an error of the size ff(e? Q s) « 10 9 
which is of no interest to us. In the case of Y/ATk, we presume the following assumptions: 



• First, we make the approximation ATk = F\ which is good up to the precision of 10 3 . 

• Second, we assume that the contribution from the two pion state is dominant in Y,T\ which 
is good in the precision level of 10~ 3 . 

• Third, we assume that the contribution from the 1 = two pion state is dominant in 
compared with that from the 1 = 2 state. This approximation is good up to the precision of 

io- 7 . 



Using these assumptions, we can approximate the Y/ATg as follows, 

Then, we can express Ek approximately as follows, 

£K = (^(/ = 0)|^|^) " £ <°) + * + " ' VAAfff ^7 A£ ^ • (2 " 10) 



The correction of ^(£ 3 j) in Eq. ( 2.10 ) is smaller than both the current experiment precision 



and the size of the long distance contributions of the m' [|7|]. The last expression in Eq. ( 2.1C ) is 



obtained by substituting £m) with Eq. ( |2.6[ ) and Eq. ( |2.9[ ). The coiTection of Ask contains both 
short-distance (SD) contribution and long-distance (LD) contribution, which are expected to be 
about « 5%[[7|]. Here, we also neglect this contribution from Ask, mainly because it is not known 
to a sufficient precision theoretically. 
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In this analysis, we take into account only the short-distance contribution from the dimension 
6 operators m',^. In the SM, this part can be calculated from the box diagram [JsJ] - Here, we follow 
the notations in [__]. Then, we can obtain the following master formula which will be used in this 
analysis: 

V2smd(c e B K X + ^o), (2.11) 



where 



i„SMi 
\ b K I 



X = fiX 2 \V,- a 



cb X 



|V c fo| 2 (l -p)r\2So(x t ) + ri 3 So(x c ,x t ) - T7iS (x c 



G 2 F F^mgoM^ 
6V2n 2 AM K ' 



(2.12) 
(2.13) 



where we use the experimental value for AMk- The input <^o has been taken from the lattice calcu- 
lation which accounts the long distance contribution. 



(a) 



G F 


1.16637(1) x 1(T 5 GeV~ 2 




1 




M w 


80.399(23) GeV 




1 






m c {m c ) 


1.25(9) GeV 




C] 


0.808(22) 




| CKMfitter 


m,(m,) 


162.7(1.3) GeV 




0.832(17) 




| UTfit 


ill 


1.43(23) 






0.2253(7) 




| CKMfitter 


m 


0.5765(65) 




C] 


0.2246(11) 




| UTfit 




0.47(4) 




C] 


u.uz_ 014 




| CKMfitter 


9 


43.51(5)° 




1 


0.130(18) 




| UTfit 


m K o 


497.614(24) MeV 




1 


0.341(13) 




| CKMfitter 


AM K 


3.483(6) x 10" 12 MeV 




] 


0.350(13) 




| UTfit 



(b) Wolfenstein Parameters 





156.1(0.2)(0.8)(0.2) MeV 








156.1(1.1) MeV 




1 1 


] LAT.AVG. 


B k 


0.7674(99) 
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] LAT.AVG. 


0.727(4)(38) 




J_SWME 


6) 


-1.63(19) (20) x 10" 4 




13] 



\Vcb\ 


41.85(42) (9) (59) x 10" 3 


I 


12 


| (X c lv+X s y) Km 


41. 68(44) (9) (58) x 10" 3 


[ 


12 


1 (X c /v) K i„ 


41.87(25) x Kr 3 


[ 


12 
12 


I (X c lv + X s y) ls 


42.31(36) x 10~ 3 


[ 


1 (X c lv) ls 


41.5(7) x 10" 3 




ljj_Incl.PDG.AVG. 


39.5(1.0) x 10~ 3 


I 


1 1 1 Excl. 



(c) Lattice Calculation 



(d) V cb 



Table 1 : Input Parameters 



3. Input Parameters 

The parameters, m c ,m t ,r]i,r\2,T\3 depend on the renormalization scheme, and so are taken 



from the single reference for consistency (Table Ha] 1 ). The CKMfitter and UTfit results in Table lb 
are obtained by their own global fit method using the same PDG inputs. 

In Bk calculation in Table BMW quotes the smallest systematic error. It dominates the 
smallness of the lattice average error. RBC-UKQCD collaboration calculates ImA2 on the lattice. 



'in Ref. |__], they reported results of T]i up to NNLO but end up with a noticeably larger error bar. Hence, we decide 
to use the NLO value. 
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Using this value, they determine <^o through the relation 



Re 



1 ReA 2 /ImA 2 \ 
VeJ ^2\e K \ ReA W 2 *V' 



Other inputs such as ReAo,ReA 2 ,£if, and e' k /£k are taken from experiments. 

Inclusive V c b can be extracted from global fit of measured moments (lepton energy, hadronic 
mass, and photon energy) of the decay channels: 

B^XJv, B^X s y. 

We use the PDG average value as the representative of the inclusive V c b- The quoted exclusive V c b 
is the average of two semi-leptonic decay channels: 

B -»• D*£v , B -> D£v . 

For each scalar and vector channel, HFAG result is combined with FNAL/MILC lattice QCD cal- 
culation of the zero recoil form factor. 



4. Error Estimate 

For the function with N arguments, the error propagation formula gives the combined error oy 
in terms of the errors of each arguments: 



T 2 £ r df(x) 

],k=\ J 



dm 

(x) dx k 



(4.1) 



where Cjt denotes the normalized correlation between the parameters Xj and x^, and \Qj\ < 1. 
Especially the diagonal components C„ = 1. We turn off the correlation and so Qj = 5/y. In the 
case of asymmetric error, p given by CKMfitter, we take a larger error and treat it as a symmetric 
error. 

For e| M , 

/(x) = |e*WI (4-2) 
x=(d,G F ,F K ,m K o,M w ,AM K ,^ ,B K ,X,p,f},\V cb \,ri l ,ri2,ri 3 ,x c ,x t ). (4.3) 



To check the unitarity of the CKM matrix, 



y; 7 (x) = [v(x)v t (x)-/] !7 = A ;7 

x = (A,X,p,ri;\V cb \). 



(4.4) 
(4.5) 



We use the Wolfenstein parametrization to evaluate each elements of the CKM matrix Vy, except 
for V c b itself. Real and imaginary part of A;,- are separately treated. 
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(a) CKMfitter (b) UTfit 



Figure 1: £a:(x10 3 ) 

5. Results 

The Wolfenstein parameter set, Bk , and V cb has a multiple choice. It forms 8 input parameter 
sets. We calculate ef M for these sets as shown in Fig. |]. We find out that |%| shows 3.1a tension 
between |£^ xp | and |£f M |, using exclusive V cb and SWME calculation of Bk as shown in Fig. [Ta|. 
With the UTfit Wolfenstein parameters, the tension is slightly reduced to 2.9a [Fig. |Tb| ]. 

The deviation matrix A provides a test of the Wolfenstein parameters, outputs of the global fit. 
These are inputs for the e K M , as well. So testing the compatibility between the global fit results and 
the V c b is needed to interpret the difference between e| M and e K M . We find that the numerical size 
of Ajj has the following hierarchy 

10V(|ReAi 2 |) = 10^(|ReA 22 |) = ^(|ReA 23 |) . (5.1) 



ReA 22 shows that the difference between exclusive V cb and Wolfenstein parameters from CKMfitter 
(UTfit) are about 1.0a (1.8a) as shown in Fig. and 2b. ReA 2 3 shows the difference about 1.0a 
(1.8a) as shown in Fig. |c| and 2d Other components of A, which does not depend on the choice 
of inclusive or exclusive V cb , are so small as to be consistent with the unitary ansatz. 

In Table ||, V cb dominates the error of e K M regardless of the inclusive or exclusive determina- 
tion. In case of the SWME calculation of Bk and inclusive V cb , both contribute to the total error in 
comparable size. In the case of the lattice average of Bk, V c b becomes an extremely dominant error 
and the subdominant error comes from f] . 



W.P. 


Vcb 


B K 


m c 


m 




F K 


B K 




P 


n 


V cb 




Incl. 


LAT.AVG. 


10.62 


5.24 


12.50 


2.36 


1.98 


1.85 


4.29 


17.29 


41.43 


CKMfitter 


SWME 


8.10 


4.00 


9.53 


1.80 


25.06 


1.57 


3.28 


13.19 


31.61 


Excl. 


LAT.AVG. 


7.94 


3.92 


9.34 


1.53 


1.28 


1.68 


2.63 


11.17 


58.99 




SWME 


6.61 


3.26 


7.78 


1.27 


17.67 


1.56 


2.19 


9.29 


49.10 



Table 2: Error Fractions of / £ ct?. The UTfit Wolfenstein parameters(W.P) show the same tendency. 
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Incl.Vcb Exc\.V C b 
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(b) ReA 22 (UTfit) 



2 
1.5 

1 
0.5 



-0.5 

-1 
-1.5 

-2 
-2.5 

-3 
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(c)ReA 23 (CKMfitter) 
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-1 
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Figure 2: ReA 22 (x 10 4 ) and ReA 23 (x 10 3 ) to test the CKM unitarity. 
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